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Introduction 

The extraordinary biconcave shape of a red blood cell has attracted much interest for 
many years. In the last two decades, people generally accepted the shape of biological 
membranes such as blood cells is closely related to the formation of lipid bilayer vesicle 
in aqueous medium. Based on the elasticity of lipid bilayers proposed by Helfrich [HI], 
the shape S, regarded as an embedded surface in M^, is determined by the minimum of 
the bending energy involving the volume V(S) enclosed by S, the area A(S), the mean 
curvature H and the Gaussian curvature of S. More precisely, Helfrich suggested to 
study the bending energy 

{2H + Co)' AA+^yA KdA + AA(S) + pV(S), 

where fee, kci Cq, A, and p are constants interpreted as follow: kc is the bending rigidity, k^. 
the Gaussian curvature modulus, Cq the spontaneous curvature, A the tensile stress, and 
p = Po — Pi the osmotic pressure difference between the outer (po) and inner (pj) media. 
Here, we have taken a geometric sign convention so that Helfrich's original bending energy 
should be written in the above form. 

According to the Gauss-Bonnet Theorem, the second integral in the bending energy 
is a topological constant. Therefore, within a certain topological class of S, it is sufficient 
to study the functional, 

^(S) = jj2H + cqY AA + AA(S) + pV(S) 
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where cq, A = 2\/kc^ and p = 2p/k^ are the constant parameters of the functional. This 
functional will be referred as the Helfrich functional in this article. 

Many interesting surfaces such as minimal surfaces, constant mean curvature surfaces 
and Willmore surfaces can be regarded as critical points of the Helfrich functional for 
suitable combinations of the parameters Cq, A, and p. New axisymmetric explicit solutions 
had also been found recently [N002] . 

When Co = A = p = 0, the functional ^(S) is referred as Willmore functional in 
differential geometry which has been widely studied in recent decades. Moreover, there 
are some important analysis and open problems concerning the Willmore functional, [W, 
ch. 7]. As it is observed by physicists [DH], Willmore functional is not a good model 
for the shape of red blood cells. This fact can also be seen by the result of geometers 
that the unique minimum of the Willmore functional for topologically spherical vesicles 
(embedded surfaces of genus zero in terms of topologist) is the round sphere. Therefore, 
not all combinations of the parameters give stationary vesicles of shape similar to red 
blood cells observed experimentally. 

On the other hand, a typical way to investigate the minimizing surface is by finding 
solutions to the variational equation of the functional J^. A large class of axisymmetric 
stationary vesicles of spherical and toroidal topology has been calculated, [DH, L, S, MB, 
OY]. In these works, among other interesting shapes of lipid bilayer vesicles, the shape of 
the red blood cell can be simulated numerically or given by a special solution a specific 
variational equation with suitable combinations of the parameters. It is natural to ask for 
the conditions on the parameters such that the Helfrich functional possesses a stationary 
vesicle of biconcave shape. 

Main Results. In this article, we are going to give clear conditions on how biconcave 
axisymmetric surfaces are formed. More precisely, we find a sufficent condition for that 
the Helfrich shape equation of axisymmetric vesicles to have solutions of biconcave shape. 
Besides, we exhibit that when the equation has a solution with reflection symmetry or 
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of biconcave shape, certain geometric quantities of the vesicle must obey some conditions 
governed by the parameters. These conditions on the parameters are very mild. In 
particular, the case that Cq > 0, A > 0, and p > is sufficient to ensure the existence of 
biconcave solution. We also briefly comment on how the combination of the parameters 
Co, A, and p affects the existence of biconcave solution. 

The sufficient condition for the formation of biconcave vesicle may be formulated in 
terms of a cubic polynomial, 

(1) Q{t)^t^ + 2cot^ + {cl + ~X)t-^. 

We proved that if all roots of Q{t) Eire positive, then there is always an ax- 
isymmetric biconcave vesicle which is the stationary surface for the Helfrich 
functional. This result will be summary by the table at the end of this section (p. 5) 
along with some typical pictures of Q{t) and the corresponding graphs of the solution and 
its derivative. 

It is easy to verify that Q{t) < — | for alH < if Cq > 0, A > 0, and p > and hence 
all roots are positive. Mathematically, that all roots of Q{t) are positive can be written 
as 

max{Q{t) : -oo < t < 0} < 0. 

Explicit formula for this in terms of the parameters can also be found. However, the 
above is more concise and precise enough. 

For the necessary condition, we observe that if Cq > and there is an axisymmetric 

stationary vesicle of biconcave shape, then its curvature at the center must be smaller 
than the first positive root of Q{t). In fact, we have 

2cow','' + {cl + X)w', - I < 0, 

where Wq is the meridinal curvature (or w'q^ is the Gaussian curvature) at the center. 

Furthermore, if the axisymmetric vesicle has a refiection symmetry with respect to a 
plane perpendicular to the rotation axis, then we obtain a relation bewteen the Gaussian 
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curvature at the "equator" (the intersection circle of the axisymmetric vesicle with the 
plane of reflection) and the radius r^o of the "equator" in terms of the polynomial Q{t) 
given by (1) as follows. 



This paper is organized as follows. In §1, we rewrite the Helfrich shape equation 
of axisymmetric vesicles into two useful forms that are more convenient for our later 
discussion. In §2, the change of the principal curvature in meridinal direction is analyzed 
and a necessary condition for the existence of biconcave surface is given. In §3, we present 
a sufficient condition for the existence of biconcave surface. In §4, we give a brief discussion 
of other situation such as axisymmetric vesicles which do not have reflection symmetry 
or which is not biconcave. 

We would like to thank K. S. Chou for his discussions and valuable suggestions on the 
writing of this manuscript. 





1. Equation for Axisymmetric Vesicles 

In this article, we will study axisymmetric solution surface E of the Helfrich variation 
problem which has a reflection symmetry by the plane perpendicular to the rotational 
axis. If the rotational axis is labelled 2;-axis and the plane of reflection xy-plane, then the 
surface E can be obtained by revolving a radial curve about the z-axis on the upper half 
plane and reflecting it to the lower half. 

0.8 
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02 
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A cross-section of a biconcave axisymmetric surface 
(with Co = 1, A = 0.25, p = I) 

In other words, one may parametrize the upper part of E by X = (r cos6', r sin6', 2;(r)), 
with a function z{r) > defined for r in some interval [0, Too], where Too is the radius 
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of the "equator" that is determined by the surface. The natural boundary conditions 
are z{0) > and z{r(X)) = 0. The rotation and reflection symmetries of the surface then 
impose the conditions that z'{0) — and z'{r) — > — oo as r — > Too- Finally, to model for a 
biconcave surface, we also need to assume that z"{0) > in order to have a solution that 
is concave at the center. Moreover, the biconcave shape of the surface also confines the 
graph of z to have a unique point of inflection. This is equivalent to require that z' has 
only a unique maximum and no other critical point. A cross-section of the upper part of 
S is shown in the picture. 



Equation and Conditions. Usually, the Helfrich shape equation of axisymmetrc vesi- 
cles is written in term of the angle ip between the surface tangent and the plane perpen- 
dicular to rotational axis [DH, JS, ZL]. However, it is convenient for our discussion to 
rewrite the equation into a equation on the derivative of the graph w{r) — z'{r) — tanip. 
Then the Helfrich shape equation of axisymmetrc vesicles becomes an equation for w, 
which is 

2r ,, 5rw ,2 2w' 
11! = 1/; 

(l + w2)5/2 (l + ^/;2)7/2 (1 + m;2)5/2 

2w + 2cqw'^ (cq + X)rw pr^ 

+ r(l + w2)3/2 + TTw2 ^ (1 + ^2)1/2 ~ ~Y ' 

We are going to study the solution w{r) to this equation with initial choice u'(O) =0 and 

^'(0) = Wq > 0. Our goal is to find solution w that also satisfies the end point condition 

that wir) —00 as r ^ Too and an integral condition that —00 7^ / wdr < 0. The 

Jo 

integral condition is equivalent to z{roo) = < z{Q) < 00 which ensures that z{r) to- 
gether with its refiection produce a closed axisymmetric surface without self-intersection. 
Furthermore, w is required to have a unique local maximum and no other critical points. 
A typical picture for the graph of such w is shown. 



(2) 
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A typically expected graph of w = z' 



There are two useful ways of grouping the lower order terms in the equation. A careful 
study of these lower order terms reveals the quahtative behaviors of the solution w. One 
of the grouping involves the polynomial Q{t) that we have mentioned in the introduction. 

Let Q be the cubic polynomial Q(t) — + 2cot^ + (cq + X)t — - and R(t) be the 
quadratic part of Q(t), that is, R{t) = Q{t) — t^. With this notation, after multiplying 
with rw', the equation (2) can be written as follow. 



(2a) 
(2b) 

where K{r) — 





/ 








/ 


_(l + «;2)5/2_ 





w 



(1 + «;2)V2 
-2 



+ r^w'R{K,{r)); or 
+ rVg(«;(r)), 



w 



is the principal curvature of E in the meridinal direction. In 



other words, the lower order terms capture the changes of the meridinal curvature in 
the differential equation. Note that if Cq = A = p = 0, the solution to equation (2a) 
corresponds to the situation that z = z{r) is a circular arc with radius 1/wq. This is 
exactly the case when the meridinal curvature is a constant. 

The initial choice w'q is actually the meridinal curvature at the center of E, i.e., 

w 



K,{r) — 



Moreover, for both Q and i?, (5(0) = i?(0) = — p/2 < 0. Therefore, if w'^ is less than the 
smallest positive root of (5(t), then both Q{wq) < and R{wq) < 0. Our analysis (in 
§3) shows that these negativity conditions are essential for the meridinal curvature «;(r) 
to change sign. We will also demonstrate that, for any set of parameters satisfying our 



0' 



as r 



0. 
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condition, a solution w{r) with Wq small enough always satisfies all of our requirements. 
In fact, if exact values of the parameters are given, one can always calculate the suitable 
value of Wq numerically. 

2. Necessary conditions for the formation of biconcave solution 

It turns out that the principal curvature of the vesicle E in the meridinal direction, 

w 

which we denoted by ^(r) = — , is the most basic quantity in our analysis. The 

rv 1 + w'^ 

equation (2b) implies that if cq > and the initial curvature at the center w'q is large, in 
the sense that R{wq) > 0, then K{r) increases and a biconcave shape cannot be formed. 
So a necessary condition for formation of axisymmetric solution of biconcave shape in the 
case that cq > is 

R{w'o) = 2cow'^^ + {cl + X)w'o - I < 0. 

The following pictures show the behaviors of K{r) for two different values of Wq with 
Co = 1, A = 0.25, and p = 1. For these parameter values, i?(0. 277124) = 0. The one on 
the left hand side is taken with Wq = 0.278 so the necessary condition is not satisfied. 
In this blows-up at finite distance from the rotational axis so the corresponding 

solution will not lead to a closed surface of the desired biconcave shape. The picture on 
the right hand side is taken with Wq = 0.276. In this case, the necessary condition is 
satisfied, but k only decays gradually and the solution w{r) still does not give a biconcave 
surface. Numerically, we see that w{r) — oo as r ^ Too but the total area is positive, 
i.e., the surface together with its refiection does not close up to form a closed surface 
in W without self-intersection. However, we can only show that R{wq) < implies the 
decrease in K,{r) without knowing whether it implies the sign change of /t(r). So, it is 
not sure whether R{w'q) < is sufficient for the formation of closed surface of biconcave 
shape, a situation in which K{r) has to change from positive to negative. 
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Graph of z{r) for a convex surface 
where R{wq) > 
Below: the corresponding graph of w = z'{r) 




Graph of z{r) leading to a surface 
with self- intersection where R{w'q) < 
Below: the corresponding graph of w = z'{r) 



0,275 
0.27 

0.265 
0.26 

0.255 
0.25 

0,245 



In order to obtained an axisymmetric solution with reflection symmetry via a solution 
of ordinary differential equation, one also needs to verify the first variation is zero with 
respect to any variation near the plane of reflection. It turns out that this is always the 
case for the solution we are going to study in the next section. Moreover, the analysis 
shows that the radius Too of the "equator" and the Guassian curvature K{roo) of the 
vesicle at any point on the "equator" must satisfy the following relation, 

-1 „ r-v 



Kiroof = -Q (-) 



This relation is easily verified for the case that cq — X — p — 0. In this case, we have 
Q{t) = and ii'(roo)^ = This is compatible with the fact that the solution surface 

^oo 

is a round sphere. 



3. Sufficient condition for the formation of biconcave solution 



In the preceding section, although that R{wq) < is known to be a necessary con- 
dition for the case that cq > 0, it may not be sufficient. Hence, we need to impose a 
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stronger condition that all roots of Q(t) are positive and require w'q to be small enough. 
Unfortunately, we do not have a general formula for the threshold value of w'q. 

The smallness assumption on w'q imphes that not only R{wq) < but also Q{wq) < 0. 

In this case, the mcridinal curvature changes sign. And in terms of w{r), though it is 

initially positive, it becomes negative when r reaches a certain value Tq. Note that for 

different initial value Wq, the value of Tq is different. Although the precise relationship is 

unknown, we know that tq is comparable to Wq^^"^ when w'q is small. More precisely, we 

have the following limiting inequalities 

16 , r2 16 
— < lim ^ < — . 

In fact, we discover that, for any point of inflection Tc of the graph of z{r) with ^(rc) > 0, 



16 



< lim 



< lim 



< 



16 



3p w'l^^O w'q w'q^Q w'q p 

As a consequence of these inequalities, the solution will have a unique point of inflection 
of the graph of z{r). Moreover, for the case that all roots of Q{t) are positive (see §2), the 
function w{r) blows down monotonically to — oo in a finite distance from the rotational 
axis after it becomes negative. 

The following are pictures of z{r) and w{r) for cq = 1, A = 0.25, p = 1, and Wq = 
0.26. This set of parameters satisfies our condition. In fact, Q{t) has a unique positive 
root at 0.268828. So, it is clear that Q{wq) < and the function w{r) goes to — oo at 
Too ~ 5.39215. However, w'q = 0.26 is not small enough to give a solution w{r) with total 
negative area and hence it does not give a vesicle of the required biconcave shape. 





A solution z which does not give an embedded surface although w = z' blows down to — oo 
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However, if we choose Wq = 0.2063860 and Wq = 0.15 for the same set of parameters 
Co = 1, A = 0.25, and p = 1, then in each case, the total area under the graph of w{r) 
is nonpositive. Hence, each of the graph of z{r) can be reflected and rotated to form a 
vesicle of biconcave shape. The following pictures are the function z{r) and w{r) for these 
values of Wq respectively. 




Graphs of z{r) giving rise to biconcave surfaces and their z' below 



The mathematical proof of this result involves more delicate study of the equation (2b) 
together with sharp estimates of Tq, Too, the rate of change of w{r) at Tq, and the integral 
of w{r). Please refer to the mathematical paper by the authors for the details of the 
analysis. 



4. Discussion 

Finally, we would like to remark on other solutions which do not correspond to bi- 
concave surfaces. Firstly, when Cq + A < 0, one may still obtain solution for nonconvex 
surface even w'q is not very small. 
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A solution which gives a nonconvex surface 



The above picture is a solution for cq = 1, A = —3, p = 2 with Wq = 2. Note that w'q is 
greater than the positive roots of R{t) and Q{t). Moreover, Q{t) has negative roots and 
thus this is beyond the situation that we discuss in previous sections. In this case, after 
w — z' becomes negative, it may not decrease all the way to negative infinity, which is 
why a second "petal" may form. 

The other type is a rotational solution which is not symmetric under reflection in 
xy-plane [DH, HDH]. The following is a solution with co = — 1,A = — 1 and p = 1. 



It can be observed that the upper part of this solution satisfies our equation (2) but 
not the area condition. Although the reflected graph is also a solution locally, it does 
not give rise to a surface without self-intersection. For this combination of parameters, 
there is a bifurcation of solution at Too to close up by an asymmetric surface. Prom our 
numerical study, it should be remarked that not all solutions in the upper quadrant may 
be closed up by an asymmetric counterpart. For example, if the upper part satisfies the 
area condition, the only possible lower part is the refiection symmetric one. This suggests 
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A solution which has no reflection symmetry 
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a possible uniqueness of biconcave solution. Furthermore, we surmise that for cq > 0, 
A > and p > 0, solution of the asymmetric type never occurs. The local minimum of 
Q{t) for i < may be the underlying obstruction of such happening. However, at the 
moment, there is not yet a mathematical proof. 
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